KKN molecule state in three-body calculation 
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A KKN system with 7=1/2 and J p = l/2 + is investigated with non-relativistic three-body 
calculations by using effective KN and KK interactions. The present investigation suggests that 
a weakly bound state for the KKN system can be formed below the two-body threshold of K 
and quasibound KN with a 40 ~ 60 MeV decay width of KKN — > nYK. This corresponds to an 
excited S baryon with J p = 1 /2 + located around 1.9 GeV. Studying the wave function of the KKN 
system obtained in this formulation, we find that the three-body bound system has a characteristic 
structure of KN+K cluster with spatial extent. 

PACS numbers: 



I. INTRODUCTION 

The study of hadron structure is one of the most impor- 
tant issues in hadron physics. Recent interest in this line 
is developed in exploring quasi-bound systems of mesons 
and baryons governed by strong interactions among the 
hadrons. One of the long-standing candidates is the 
A(1405) resonance considered as a quasi-bound state of 
KN system [l[. It has been also suggested that the 
/o(980) scalar meson is a molecular state of KK Q. For 
nuclear systems, bound states of an eta meson in nu- 
clei and an anti-kaon in light nuclei were predicted in 
Refs. @, H @. Recently a multi-hadron state was pro- 
posed for the + baryon resonance to explain its narrow 
width [fj. 

In such multi-hadron systems, anti-kaon plays a unique 
and important role due to its heavy mass and Nambu- 
Goldstone boson nature. The heavier kaon mass indi- 
cates stronger s-wave interactions around the threshold 
than those for pion according to the chiral effective the- 
ory. In addition, since typical kaon kinetic energy in 
the system estimated by range of hadronic interaction is 
small in comparison with the kaon mass, we may treat 
kaons in multi-hadron systems in many-body formula- 
tions. 

The strong attraction in the KN system led to the 
idea of deeply bound kaonic states in light nuclei, such 
as K ~pp and K~ppn, by Akaishi and Yamazaki 0, 0, 
H, H E}. Later, many theoretical studies on the struc- 
ture of the K~~pp system have been done, for example in 
Faddeev calculations [ll], [13] and in variational calcula- 
tions [HI, EH, having turned out that the K~pp system is 
bound with a large width. Ex per imental search for these 
states has been reported (IH Il6l. [i~7l [l8| . whereas clear 
evidences are not obtained yet and interpretations of the 
experimental data are controversial 0, H3] . Motivated 
by the strong KN{I = 0) attractions, quests of vm\\t\-K 
nuclei are challenging issues @, HU, [H, H, Hj] . This is 
also one of the key subjects related to kaon condensation 
in dense nuclear matter [H, H3, HH • 

The key issue for the study of the KN interaction is 
the resonance position of the A(1405) in the KN scat- 
tering amplitude. The A(1405) is observed around 1405 



MeV in the 7rE final state interaction, as summarized in 
the particle data group [29j. Based on this fact, a phe- 
nomenological effective KN potential (AY potential) was 
derived in Refs. [H, [l3| , having relatively strong attrac- 
tion in the 1 = channel to provide the K~p bound state 
at 1405 MeV. Recent theoretical studies of the A(1405) 
based on chiral unitary approach have indicated that the 
A(1405) is described as a superposition of two pole states 
and one of the state is considered to be a KN quasibound 
state embedded in the strongly interacting 7rX contin- 
uum [Mil 12, 11, 11 11]. This double-pole conjecture 
suggests that the resonance position in the KN scatter- 
ing amplitudes with I = is around 1420 MeV, which is 
higher than the energy position of the nominal A(1405) 
resonance. Based on this chiral SU(3) coupled-channel 
dynamics, Hyodo and Weise have derived another effec- 
tive KN potential_(HW potential) [II]. The HW po- 
tential provides a KN quasibound state at ~ 1420 MeV 
instead of 1405 MeV, and it is not as strong as the AY 
potential. 

The common achievement of these studies on the s- 
wave KN effective potential is that the KN interac- 
tion with I = is strongly attractive and describes the 
A(1405) resonance as a KN quasibound state. Basing 
on these strong KN interactions, we examine possible 
bound states of a lightest two-anti-kaon nuclear system, 
namely KKN with 1=1/2 and J p = 1/2+, in the 
hadronic molecule picture. We expect that the three- 
body KKN system will form a bound state due to the 
strong attraction in the two-body KN subsystem pro- 
ducing the 10 ~ 30 MeV binding energy. The questions 
raised here are whether the strong attraction is enough 
for a three-body bound state below the threshold of the 
K and the KN quasibound state, and what structure the 
bound state has, if it is formed. 

The KKN molecule state may have characteristic de- 
cay patterns depending on the binding energy. If the 
three-body state is above the threshold of the K and the 
KN quasibound state, it can decay into K and A(1405), 
and, consequently, the bound state has a large width. 
If the KKN system is bound below the threshold, the 
bound state has a comparable width to the A(1405), 
and the main decay mode is a three-body ttT,K and 
two-hadron decays are strongly suppressed in contrast 



2 



to usual excited baryons. For deeply bound KKN sys- 
tems, the molecular picture may be broken down and 
two-hadron decays may be favored. Presently experi- 
mental data for S = — 2 channel are very poor. But we 
expect that such an interesting molecule state could be 
observed as an exited state of S baryon, for example, 
in double strangeness exchange (K~.K + ) processes at 
J-PARC. 

In this paper, we investigate the KKN system with 
1=1/2 and J p = 1/2+ with non-relativistic three-body 
calculations by using the HW and AY potential as effec- 
tive KN interactions. In Sec. HU we describe the frame- 
work of the present calculations. We apply a variational 
approach with a Gaussian expansion method [37j to solve 
the Schrodinger equation of the three-body system. By 
treating the imaginary potentials perturbatively, we find 
the KKN quasibound state. In Sec. IIII| we present our 
results of the three-body calculation. In analysis of the 
wave functions, we discuss the structure and the binding 
mechanism of the KKN state. The effects of KK in- 
teractions on the KKN system are also discussed. Sec- 
tion [IV] is devoted to summary and concluding remarks. 



II. FRAMEWORK 

We consider a non-relativistic three-body potential 
model for the KKN system. We calculate the KKN 
wave function by using effective two-body interac- 
tions [H, [36| in local potential forms. We apply a varia- 
tional approach with a Gaussian expansion method [37j 
in solving the Schrodinger equation for the three-body 
system. In this section we explain the details of the 
framework and interactions used in this work. 



A. Hamiltonian 

In the present formulation, we use the Hamiltonian for 
the KKN system given by 



H 



: T + V kN {n) + V RN {r 2 ) + V RR (r 3 ), 



(1) 



with the kinetic energy T, the effective KN interaction 
V Rn and the KK interaction V R k- These interactions 
are given by local potentials as functions of K-N dis- 
tances, ri, r 2 , and the K-K distance, which are de- 
fined by r x = |x 2 -x 3 |, r 2 = |x 3 -Xi| and r 3 = |x2-x x |. 
The vectors Xj., X2, X3 denote the spatial coordinates of 
the first anti-kaon (-Ki), the second anti-kaon (K2) and 
the nucleon, respectively. For convenience, we use Ja- 
cobian coordinates, r c and R c , in three rearrangement 
channels c = 1,2,3 as shown in Fig. [TJ In the present 
calculation, we assume isospin symmetry in the effective 
interactions, and we also neglect the mass differences be- 
tween K~ and K° and between proton and neutron by 
using the averaged masses. We do not consider three- 
body forces nor transitions to two-hadron decays. 
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FIG. 1: Three Jacobian coordinates of the KKN system. 



The kinetic energy T is simply expressed by the Jaco- 
bian coordinates for a rearrangement channel as 



T 



-1 



-1 



2/i r 



V 2 



(2) 



with the reduced masses (i ra and [1r c for the correspond- 
ing configuration. For instance, \i ri = MkMn /{Mk + 
M N ) and [i Rl = M K (M K + M N )/(2M K + M N ) for the 
rearrangement channel c = 1. Here and Mx denote 
the averaged nucleon and kaon masses, respectively, as 
M N = 938.9 MeV and M K = 495.7 MeV. 

For the effective interactions, V RN and V RR , we use 
local potentials obtained by s-wave two-body scattering 
amplitudes with isospin symmetry. In the present three- 
body calculations, we take /-independent potentials for 
simplicity. The details of the local potentials will be dis- 
cussed in Sec. Ill Bl The coupled-channel effects of the 
KN to other relevant channels have been already im- 
plemented to the effective single-channel KN interaction 
Vrn- Consequently V RN has an imaginary part owing 
to scattering states below the KN threshold, and the 
present Hamiltonian is not hermitian. In the calcula- 
tions of KKN wave functions, we first use only the real 
part of V RN in a variational approach, and then we cal- 
culate the energy E with the expectation value of the 
total Hamiltonian JlJ with respect to the obtained wave 
functions. The width of the bound state is evaluated by 
the imaginary part of the energy as T = — ZhaE. The 
details of the calculational procedure will be described 
later. 



B. Effective interactions 

In this subsection, we explain the details of the ef- 
fective interactions of the KN and KK two-body sub- 
systems in our formulation. We use two effective KN 
interactions which were derived in different ways, and 
compare the results obtained by these interactions to es- 
timate theoretical uncertainties. Both KN interactions 
have so strong attraction as to provide the A(1405) as 
a quasibound state of the KN system. The important 
difference is the binding energy of the KN system, as 
already mentioned in introduction. 
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One of the KN interactions which we use here is given 
by Hyodo and Weise in Ref. [36]]. This effective interac- 
tion was derived based on the chiral unitary approach [38[ 
for s-wave scattering amplitude with strangeness S = — 1 
of the lowest-lying octet meson and baryon. A equivalent 
single-channel KN amplitude can be obtained by reduc- 
ing the four (five)-channel problem for I = (I = 1) to 
a KN single-channel problem including dynamics of the 
rest channels. The local KN potential was constructed 
in coordinate space so as to reproduce the single-channel 
KN interaction as a solution of the Schrodinger equation 
for the KN system with the local potential. 

The potential is written in a one-range Gaussian form 

as 



*KN 



U RN ex P 



, , . l ir/b) 2 ]P RN (I = Q) 
+t/i= 1 exp [-{r/bf} P RN (I = 1), 



(3) 



with the isospin projection operator P R jq{I — 0, 1) and 
the range parameter b — 0.47 fm. The range param- 
eter was optimized for the parametrization referred as 
HNJH [39] in Ref. [3f|. We refer this potential as "HW- 
HNJH potential". The strength U Rn (I = 0,1) has en- 
ergy dependence and is parametrized in terms of a third 
order polynomial in the energy uj: 



Tr I=0,l/ \ 



K 



1=0,1 



„ 7=0,1 

K 1 uj - 



„j=0,l 2 

K 2 u 



K 3 UJ , 



(4) 



for 1300 MeV < uj < 1450 MeV. The coefficients K\ 
for 1 = and 1=1 are given in Table IV and V of 
Ref. (36| , respectively. Here we use the corrected version 
of the local potentials [36| . 

In chiral unitary approaches for the meson-baryon in- 
teractions, two poles are generated at 1.4 GeV region, 
and the A(1405) resonance is described as a KN qua- 
sibound state in the strongly interacting 7rS contin- 
uum [lO, In this case, the peak position of the KN 
scattering amplitudes in the 1 = channel appears at 
uj ~ 1420 MeV, which is higher than the energy posi- 
tion of the A(1405) resonance observed in the final state 
interaction of the 7r£ channel. The HW potential re- 
produces this feature of KN scattering amplitudes, and 
therefore the position of the KN quasibound state is lo- 
cated at uj ~ 1420 MeV. With the HW-HNJH potential, 
A(1405) can be calculated as a KN quasibound state by 
solving the KN two-body Schrodinger equation in the 
s-wave 1 = channel. We get a KN quasibound state 
at 1423 MeV with a width T = 44 MeV, treating the 
imaginary part of the potentials perturbatively. This is 
a weakly bound KN state located at 11 MeV below the 
KN threshold (M K +M N = 1434 MeV). The root-mean- 
square of the K-N distance (d RN ) is also calculated as 
1.9 fm. 

The HW-HNJH potential V^'^uj) depends on the 
energy of the subsystem KN. We regard, however, uj as 
an interaction parameter and use fixed values of uj, since 
the energy dependence is small in the region of interest, 



uj = 1400 ~ 1440 MeV. We take two values of u = M N + 
M K - Suj with Suj =0_and 11 MeV. The choice of 5uj = 
MeV corresponds to KN at the threshold, while Suj = 11 
MeV is for the binding energy B(KN) = 11 MeV of the 
KN bound state obtained by the HW-HNJH potential. 
In each of the fixed HW potential, we calculate the energy 
E Rc = —B(KKN) of the KKN system and also the 
energy —B(KN) of the KN system. 

The other KN interaction which we use is given 
by Akaishi and Yamazaki (AY) in Refs. [!, [Uj]. The 
AY potential was derived in a phenomcnological way 
to start with the ansatz that the A(1405) resonance 
is a K~p bound state at 1405 MeV as reported by 
PDG |29| . The AY potential is energy independent and 
was parametrized so as to reproduce the KN quasibound 
state at the PDG values of the A(1405): 

V£=°(r) = (-595 - 83i) exp [-(r/0.66fm) 2 ] , (5) 
Vj=^{r) = (-175- 105i)exp [-(r/0.66fm) 2 ] . (6) 

The AY potential has a stronger attraction in the 1 = 
channel because it was adjusted to generate a KN(I = 0) 
quasibound state at ~ 1405 MeV. It corresponds to the 
KN binding energy B(KN) ~ 30 MeV which is much 
larger than B(KN) = 11 MeV for the HW-HNJH poten- 
tial. By solving the KN two-body Schrodinger with a 
perturbative treatment of the imaginary potential, we 
get a KN (I = 0) state at 1403 MeV with a width 
of 40 MeV, and the root-mean-square K-N distance is 
d RN = 1.4 fm. Because of the strong attraction in the 
1 = channel of the AY potential, the KN{I = 0) qua- 
sibound state has the deeper binding and smaller radius 
than those given by the HW-HNJH potential. We com- 
ment that the size of the K-N state is sensitive to the 
binding energy [401 ]. 

For the KK interactions, there are few experimental 
data. The low energy theorem based on the current alge- 
bra suggests a repulsive interaction. For the s-wave inter- 
action of KK, the I = is forbidden due to Einstein-Bosc 
statistics. Thus we assume V^j^ = 0. We introduce the 
effective interaction of the subsystem KK with 1=1, 
Vkj£, in a Gaussian form 

VRR(T) = U R =i exp [~{r/b) 2 } P RR (I = 1), (7) 

where the range parameter b is chosen to be the same 
value as that of the KN interaction. The interaction 
^Rii ^ s a rea ^ f unc tion, since there are no decay channels 
open for the KK system. 

The strength U RR is estimated by theoretical calcula- 
tions of the scattering length of K + K + , which is equiv- 
alent to that of KK with 1=1. Recently, the K+K+ 
scattering length has been obtained in lattice QCD cal- 
culation as a K+K + = -0.141 ± 0.006 fm [Uj]. This value 
is consistent with the leading order calculation of the chi- 
ral perturbation theory, o^+}f+ = —0.147 fm, which is 
obtained by 



a K + K -t 



m K 

167T/I' 



(8) 
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with Jk — 115 MeV. In the present calculation, the body systems given in Ref. [37|. The spatial wave func- 
tions $j (r c , R c ) of the subcomponent of Eq. (p~0|) are 



strength U R ~ R is adjusted to reproduce the scattering 
length a££ = -0.14 fm. For the HW-HNJH poten- 
tial, with the interaction range b = 0.47 fm, we obtain 
Ugg = 313 MeV, and we find U l £± = 104 MeV with 
b = 0.66 fm for the AY potential. We will also try a 
weaker repulsion with a T £ R = —0.10 using Uhg = 205 
MeV and 70 MeV for the HW-HNJH and the AY po- 
tentials, respectively. In order to examine how strong 
the KN interaction for the three body system, at first 
we neglect the KK repulsion, then we see the effect of 
the repulsive interaction with including the effective KK 
interaction given above. 



C. Three-body wave function 

The three-body KKN wave function ^ is described as 
a linear combination of amplitudes $j (r c , R c ) of two 
rearrangement channels c = 1,3 (Fig. [ljj, since the chan- 
nel c = 2 is included by symmetrization of the two anti- 
kaon in the wave function. In the present calculation, we 
take the model space limited to l c = and L c = of the 
orbital-angular momenta for the Jacobian coordinates r c 
and R c in the channel c. Then the wave function of the 
KKN system with 7=1/2 and J p = l/2 + is written as 



l + Pv 



-4. 



V2 

*%?=o(*iM [[KK] lKK=0 N] /=1/2 
^UruRi)[[KK] lKK=1 N] I=1/2 
C=i( r 3^ R 3) [\KK] lKK=1 N] I=1 



(9) 



/2 



(10) 



where P\% is the exchange operator between the 
two anti-kaons, K\ and K 2 for two bosons. The 
\[KK]i KK N~\ I=1 j 2 specifies the isospin configuration of 

the wave function $^ (r c ,R c ), meaning that the to- 
tal isospin 1/2 for the KKN system is given by combi- 
nation of total isospin I RR for the KK subsystem and 
isospin 1/2 for the nucleon. The isospin configuration 
[[KK]i KK —oN\ / _ 1 ^ 9 @in the c = 3 is not necessary, be- 
cause it vanishes after the symmetrization in the case of 

h = 0. 

As mentioned above, we omit basis wave functions with 
l c > 1 and L c > 1 in each rearrangement channel. This 
is consistent with the fact that the effective local po- 
tentials used in the present calculations are derived to 
reproduce the s-wave two-body dynamics. We comment 
that components with non-zero angular momenta of two- 
body subsystems are contained in the model wave func- 
tion through rearrangement of three-body configurations, 
although the s-wave component is expected to be domi- 
nant. 

In solving the Schrodinger equation for the KKN sys- 
tem, we adopt the Gaussian expansion method for three- 



expanded in terms of the Gaussian basis functions, <f>^ (r) 
and ^(R): 



(r c ,R c 



max,JV m , 

E 

n c ,N C 



A 



cJkk 



,(r c )^(R c ). (11) 



The coefficients A^ 1 ^ are determined by variational 
principle when we solve the Schrodinger equation. In 
Eq. cm, and iV max are the numbers of the Gaus- 
sian basis, and the basis functions are defined by 



V^(R) = Af N e 



-\ N R A 



(12) 
(13) 



where the normalization constants are given by J\f n = 
2(2;/„) 3/4 7r- 1 / 4 and Af N = 2(2A A r) 3 / 4 7r- 1 /4 , and the 
Gaussian ranges, v n and Aat, are given by 



' n — ' mm 



Ajv = 1/-Rjv> = Rn 



Rr, 



(14) 
(15) 



We take enough bases for the present system by using 
the values given in Table U for the basis numbers and 
range parameters in the channel c = 1 and c = 3. We 
find that the mixing effect of the rearrangement channel 
c = 3 is very small in the present results of the KKN 
system. This is because the KK interactions are not 
attractive and therefore the K-K correlation is not strong 
in the KKN system. This is different from the case of 
the K~pp system where the p-p correlation is significant 
because of the attractive NN interaction. 



TABLE I: The numbers and range parameters of the basis 
functions for the rearrangement channels, c = 1 and c = 3. 



channel 


7^ max 


^"min 


'"max 


iV max 


-Rmin 


^raax 






(fm) 


(fm) 




(fm) 


(fm) 


c = 1 


15 


0.2 


20 


25 


0.2 


200 


c = 3 


15 


0.2 


20 


15 


0.2 


20 



D. Procedure of calculations 

The wave function of the KKN system is obtained by 
solving the Schrodinger equation: 

[T + V RN (n) + V RN (r 2 ) + V RR (r 3 ) - E] * = 0. (16) 

The effective interaction V RN is complex due to the pres- 
ence of the decay channels below the threshold, while 
V RR is expressed by real numbers. 
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In order to solve this equation with variational princi- of the Hamiltonian, H Rc , with variational principle in 



pie, we treat the imaginary part of the potentials pertur- 
batively. Separating the real part of the Hamiltonian, we 
write 

H Rc =T + Re Wi) + ^ W 2 ) + V RR (r 3 ). (17) 
We first calculate the wave function for the real part 
I 



the model space of the Gaussian expansion described in 
Sec. Ill CI This is equivalent to determine the eigenenergy 



E Rc and the coefficients A[ 



of the Gaussian wave 



functions <j>° c (r c )tp% c (r c ) [[KK) Ikk iV] J=1/2 by diagonal- 
izing the norm matrix and Hamiltonian matrix 



,W4(r c ) [[KK]r KK N] \H**\ 0° (r c )^ e (r c ) [[KK] Ik[ 



N 



1=1/2 



(18) 



After this variational calculation, we take the lowest- 
energy solution for H R °. The binding energy B(KKN) 
of the three-body system is given as B(KKN) — — i? Re . 
It should be checked if the bound state of the KKN is 
lower than the threshold of the K and the quasibound 
state of KN, because the KKN solution obtained above 
the KN + K threshold is not a three-body bound state 
but a two-body continuum state with K and KN quasi- 
bound state. 

Next we estimate the imaginary part of the energy E 
for the total Hamiltonian H by calculating the expec- 
tation value with the wave function W obtained by the 
Hamiltonian H Rc : 



= (V\lmV RN \V) 



(19) 



The total decay width for KKN is estimated as T — 
—2E lm . In the present calculation, we have only three- 
body decays to ttT,K and ttAK by the model setting. 

We also calculate several quantities characterizing the 
structure of the three-body system, such as spatial con- 
figurations of the constituent particles and probabilities 
to have specific isospin configurations. These values are 
calculated as expectation values of the wave functions. 
The root-mean-square (r.m.s.) radius of the K distribu- 
tion is defined as the average of the distribution of each 
anti-kaon by 



^J^IKxf+xDI*), 



which is measured from the center of mass of the three- 
body system. For the two-body KN system, r K is given 
by the spatial coordinate of K, x K , measured from the 



center of mass of the two-body system as r K — (x^.). 

We also calculate the r.m.s. value of the relative K-K 
distance defined by 



d-KK 



(*|r§|tt) 



(21) 



The probabilities for the three-body system to have the 
isospin Irk states are introduced by 



n([ 



KK 



= {^\Pkk{IkkW 



(22) 



where Prr(^kk) is the projection operator for the 



isospin configuration [[-fT-K] 



-N 



, as given before. 



' J 7 = i/ 2 ' 

We also calculate the r.m.s. radius of the K distribution 
for the each isospin state 



\1$\\{A+A)\Pkk{Ikk)*) 



{*\Prk{Ikk)\*) 



(23) 



which is normalized by Eq. (|22|) . 

In order to investigate the structure of the KKN sys- 
tem further, we calculate the expectation values with the 
unsymmetrized wave function $ given in Eq. (|10|) . Al- 
though these expectation values are not observable in 
the real KKN system having two identical bosons, they 
are helpful to understand the structure of the three-body 
system and to investigate the symmetrization effect. Wc 
calculate norms of the wave functions for the isospin 
\[K2N]i KN Ki\ iI _ 1 ^ 2 states, where the total isospin 1/2 

for the KKN system is given by combination of total 
isospin Ik n for the K^N subsystem and isospin 1/2 for 
the Ki, 



n 



(<S>\Pk 2 n(Ikn)\$) 



(24) 



where Pk 2 n(^kn) is again the isospin projection opera- 
tor. Note that this cannot be interpreted as a probability 
(20) for the [[K 2 N] Ikn K 1 ] /=1/2 states in the physical KKN 

system, since the norm is not normalized to be unity 
without the symmetrization of the anti-kaon wave func- 
tions. The r.m.s. values of the K 2 -N distance and the 
K 2 N-K\ distance for $ are evaluated as 



l (K 2 N)-K 1 




(25a) 
(25b) 



respectively. 

The perturbative treatment performed above is jus- 
tified qualitatively in the case of |(\P|ImV|\I')| <C 
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|(*|Rey|*)|. For the two-body system KN, we 
find the perturbative treatment good, observing that 
\(lmV RN ) \ = 22 MeV is much smaller than \{ReVjc N )\ ~ 
100 MeV in the HW-HNJH potential case, for instance. 
This is responsible for that we get the reasonable energy 
E = 1423 - 22i MeV for the KN with the perturba- 
tive treatment. (In a full calculation with the KN effec- 
tive interaction, the scattering amplitude reproduces the 
A(1405) resonance_at 1420 MeV with 40 MeV width) Also 
in the case of the KKN system, it is found that the abso- 
lute values of the perturbative energy | 1 Im V/f jv 1^)1 = 
20 ~ 30 MeV is much smaller than the real potential en- 



ergy |<*|ReVj 



KN 



ReV KK \V)\ = 100 



200 MeV in the 



present calculations. This is because the dominant con- 
tribution of (Vjj N ) comes from 1 = channel which has 
the strong attractive potential with the weak imaginary 
part compared with the real part. 



III. RESULTS 

In this section, we show the results of calculation of 
the KKN system with 7=1/2 and J p = 1/2+. For the 
KN interactions, we compare two potential, the HW- 
HNJH potential and the AY potential, as discussed in the 
previous sections. We first show the calculation without 
the KK interactions VkKi m or der to see if the KN 
interaction is strong enough for binding the three-body 
system. Later in subsection IIII B| we discuss the effect of 
Vrk by introducing possible repulsive KK interactions 
in the 1=1 channel. 



A. Bound KKN state without KK interaction 

In this subsection, we present the results of the KKN 
state calculated without the Vkk interaction. 



larger than that of the two-body KN bound state. It in- 
dicates that the KKN state is loosely bound with a large 
radius. In the case of the AY potential, the KKN energy 
is -36 MeV, which is 5 MeV below the KN+K thresh- 
old (—31 MeV). The AY potential gives a deeper binding 
and a smaller radius, r K = 2 fm than the HW-HNJH 
potential, reflecting the stronger KN (I = 0) attraction 
in the AY potential. 

The KKN -> ttYK decay width of the KKN_ state is 
evaluated by the imaginary KN potentials as T(KKN — > 
ttYK) = — 2E lm . In the present results, we obtain the 
width in the range of 51 MeV to 57 MeV. This implies 
that the KKN state has a comparable width to that 
of the KN state for A(1405). It is interesting to note 
that the dominant contribution in the decay width of the 
KKN state comes from the 1 = channel for the KN 
subsystem, while the 1=1 channel gives only a few MeV 
contribution. This means that the KKN state has domi- 
nant -kYiK decay and relatively small ttKK decay modes. 
This characteristic decay pattern comes from the strong 
attraction in the KN channel with 1 = reproducing 
the A(1405) as a quasibound state. 



K+K+N threshold 



> 



-10 



£Q -20 

a 
W 

-30 



-40 



(a)HW 



.At405+#- 

KKN 



(b)AY ; 



1. Energy, width and decay modes of KKN state 

First of all, it is very interesting that the KKN bound 
state is obtained below the threshold of the K and the 
quasibound KN state in both calculations with the HW- 
HNJH and with the AY potentials, as seen in Fig. [21 
where we show the level structure of the KKN sys- 
tem. In Table HI we show the results for the energies 
and radius of the KKN state as well as those for the 
KN state. For the energy-dependent HW-HNJH po- 
tential, we take two energies 5u> = and 11 MeV with 
lu = M n + M k -Suj. In the table, B(KN) and B(KKN) 
denote the binding energies for the KN and KKN sys- 
tems measured by the energies of the two-body and three- 
body break-up states, respectively. 

The KKN bound state appears as small as 1 MeV 
below the threshold of the K and quasibound KN state 
(-B{KN)) in both cases of the HW-HNJH potential. 
The calculated value for is about 4 fm, which is much 



FIG. 2: Energies of the KKN calculated with (a) the HW- 
HNJH (5ui = 11 MeV) potential and (b) the AY potential 
without the KK interaction. The calculated thresholds of K 
and the quasibound KN state are also shown. The K+K+N 
threshold is located at 1930 MeV. 



2. Structure of KKN state 

Let us discuss the binding mechanism and the struc- 
ture of the KKN system with / = 1/2 in detail. For this 
purpose, we analyze the wave functions obtained in the 
present few-body calculation in terms of the isospin con- 
figuration and the spatial structure of the KKN system, 
since the potential energy is determined by both of the 
isospin and spatial structure. 

We hrst discuss the binding mechanism based on the 
isospin configurations. For convenience, we introduce the 
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following notations; [KK]j denotes the isospin configu- 
ration \[KK] Ikk=i N'\ i _ 1 ^ 2 for the two-body KK sub- 
system with I — or 1 in the three-body system, and 
[KN] stands for [[KN] Ikn=q K] /=1/2 for the two-body 

KN subsystem with 1 = 0, which is A(1405) + AT-like 
isospin configuration. 

According to the following argument, it is easy to un- 
derstand that the [-K"-K"]i and [KN]o configurations are 
energetically favoured to gain potential energy of the 
KN subsystem. The effective KN potential Vk N has 
strong attraction in the I — and weak attraction in the 
1=1 channel. From group-theoretical arguments it is 
easily found that, on one hand, the former configuration 
[KA']i consists of [KN]q and [KN]i components with a 
ratio 3:1. As a result, both of two anti-kaons effectively 
feel the potential jVj~^ + \Vk^, which is moderate at- 
traction. On the other hand, in the latter configuration 
[KN]o, one of the anti-kaons couples with the nucleon in 
Ikn = and gains the strong attraction of V^t9, and, 
at the same time, the other anti-kaon feels much weaker 
attraction as \Vj=® + fV^ 1 . 

In the present calculation, it is found that the prob- 
ability n([AT!T]i) for the three-body system to have the 
[KK] i configuration is dominant in the KKN wave func- 
tions as shown in Table HTT1 in which II (\KK]i) is 0.87 
for the result of the HW-HNJH (Suj = 11 MeV) poten- 
tial, and it is 0.91 in the case of the AY potential. These 
values are in between two limits, [KN]q and [-K\K]i; In 
the [KN] limit, the probability II([M]i) should be 
0.75 when the symmetrization of two anti-kaons is ig- 
nored, while II([_R'A']i) = 1 for the pure \KK]i state. 
The present calculation implies that the KKN state is 
regarded as an admixture of the isospin configurations 
[KK]i and [KN)o. Therefore, the rearrangement of the 
isospin configurations is essential in the KKN bound 
state. 

Next we discuss the spatial structure of the KKN 
state. Since two anti-kaons are identical bosonic parti- 
cles, K\ and K2 cannot be identified in the symmetrized 
wave function For an intuitive understanding, it is 
helpful to analyze the wave function <I> obtained before 
the symmetrization. For the wave function <£>, we can 
definitely calculate the expectation values for the K2-N 
and (K2N)-Ki distances in the rearrangement channel 
c = 1 as given by Eq. (|25p. In the calculated results, 
it is found that N is almost the same as the K-N 
distance in the KN quasibound state, while d^ 2 N)-Ki 
is remarkably large. Actually, d(_K- 2 -/v)--R"i is more than 
three times larger than d^ ^ N i n all choices of the inter- 
actions as shown in Table. ITTTI For example, we obtain 
d-K 2 N — 1-6 f m and d(_ff 2 Af)-i?i = 6.2 fm with the HW- 
HNJH [5lo=11 MeV) potential. It indicates that one 
of the kaons widely distributes around the nucleon with 
very loosely binding and the other kaon is moving in the 
vicinity of the nucleon. In addition, the wave function 
$ contains dominantly the [A^A^Jo component and little 
the [K2N]i configuration, which are shown as the norms 



of the corresponding wave functions in Tabl dllll This 
means that the K% N subsystem has more likely the 1 = 
component and, thus the <E> has a A(1405) + K\ cluster 
structure [lij. It is worth noting that, even though the 
group theoretical argument suggests that the probability 
n([A'A^]i) for the normalized wave function should be 
0.75 without the symmetrization, the [A'A'Ji component 
is obtained as n([ATAT]i) ~ 0.9 for the symmetrized wave 
function 'J, as discussed above. This is because of the 
symmetrization effect. 

Now let us combine the analysis of the isospin and 
spatial configurations. The KKN bound state can be 
interpreted as a hybrid of the two configurations: the 
[KK]i in the inner region (I) and the A(1405) + A' cluster 
in the asymptotic region (II) as shown in the schematic 
figure (Fig. [3]). As mentioned above, one of the anti- 
kaons distributes in a spatially wide region and the other 
anti-kaon distributes near the nucleon. In the inner re- 
gion (I), two anti-kaons are coupled to isospin symmetric 
with the isospin configuration [KK]i because they oc- 
cupy the same orbit and are spatially symmetric. In the 
outer region (II), the nucleon and one of the anti-kaons 
form a A(1405) state and the other anti-kaon is mov- 
ing around the A(1405). This kaon-halo like structure 
is reminiscent of the neutron-halo observed in unstable 
nuclei. 



B. effect of repulsive KK interaction 

In the previous subsection, we presented the results 
calculated without KK interactions. We discuss pos- 
sible effects of the KK interactions. As already men- 
tioned, the I = is forbidden for s-wave KK states due 
to Einstein-Bose statistics. We here investigate how the 
repulsive KK interactions may affect the KKN state by 
introducing for the 1=1 channel. 

We use the Vfff> interactions adjusted to reproduce 
the scattering length a 1 ^ = —0.14 fm of lattice QCD 
calculation. We also use a weaker repulsion with = 
—0.10 fm. The calculated results with Vj^ are shown 
in Table [TVj. It is found that the KKN bound state 
is obtained even with the possible repulsive KK inter- 
actions. Because of the repulsive KK interactions, the 
A(1405) + K cluster develops more and the anti-kaon is 
further loosely bound. As a result, the absolute value 
of the imaginary energy lmE decreases as shown in Ta- 
ble |IV] The value of d^ 2N ^_g; 1 is extremely large as 
d(K 2 N)-R 1 — 16 fm in the HW-HNJH potential, and the 
A(1405) + K cluster feature becomes more remarkable. 

There still remains ambiguity of the strengths of the 
KK interactions due to few experimental data. Here we 
estimate possible boundaries of the repulsive interaction 
for the formation of the bound KKN state. For the 
AY potential, the KKN state is still bound even with a 
further strong KK interaction as ajfj£ = —0.20 fm. In 
the case of the HW-HNJH potential, the KKN energy is 
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TABLE II: Energies and root-mean-square (r.m.s.) radii and distances of the KN and the KKN states calculated without 
the KK interaction. For the KN interactions, the AY potential and the HW-HNJH potential with the energy parameter 
u> = Mjv + Mk — 5u> (MeV) are used. For the KN state, the real energy E — —B(KN), the imaginary energy E Irn , the 
r.m.s. K-N distances(djf JV ), rn are shown. For the KKN state, the real energy E = — B(KKN), the energy relative to the 
KN + K threshold, the imaginary energy E Iln , are shown. We also list the expectation values of ImVj-^ and ImV^ 1 
separately. 





AY 


HW-HNJH 








5io = 


Slu 


= 11 


KN{I = 


0) state 








-B(KN) (MeV) 


-30.6 


-10.4 




-11.4 


E lm (MeV) 


-19.9 


-20.9 




-21.8 


dRN (fm) 


1.4 


2.0 




1.9 


r R (fm) 


0.9 


1.3 




1.2 


KKN(I = 


1/2) state 








-B(KKN) (MeV) 


-36.0 


-11.3 




-12.6 


—B(KKN) + B(KN) (MeV) -5.4 


-1.0 




-1.2 


E lm (MeV) 


-28.3 


-25.4 




-26.8 


(ImV I=0 ) (MeV) 


-24.8 


-24.2 




-25.5 


(ImV /=1 ) (MeV) 


-3.5 


-1.2 




-1.3 


r R (fm) 


2.0 


4.2 




3.8 



TABLE III: Properties such as the isospin configurations and radii in the three-body KKN system calculated without the 
KK interaction. In the upper part, we show the values in the total wave function $ obtained after the symmetrization. The 
probabilities P([K K]o,i) to have the [.K\K']o,i configurations are listed. The r R in the [KK]o and [-^-K]i components are shown 
separately. The calculated values for d R K are also shown. We also present the expectation values for the unsymmetrized 
wave function $. The norms of the wave functions for the isospin configurations [A^VJo and [KzN]i, and distances d R2N and 
d(R 2 N)-K 1 calculated with the wave function $ are shown. 





AY 


HW-HNJH 






Slu = 




= 11 


expectation values for * 


T1([KK} ) 


0.09 


0.13 




0.13 


RdKK],) 


0.91 


0.87 




0.87 


t k\ikk=o (fm) 


2.7 


5.2 




4.8 


t k\ikk = i (fm) 


1.9 


4.0 




3.7 


dan (fm) 


3.1 


6.4 




5.9 


expectation values for $ 


n([^ 2 iV]o)$ 


0.80 


0.81 




0.80 


TL([K a N]i)* 


0.00 


0.00 




0.00 


dR 2 N (fm) 


1.2 


1.7 




1.6 


d(R 2 N)-K 1 (fm) 


3.2 


6.8 




6.2 



found above the KN + K threshold with more repulsive 
KK interactions than a 1 ^^- = —0.15 fm. In such a case, 
a resonance state or a virtual state may appear near the 
threshold. This might have a larger width than that of 
the three-body bound state. For the detailed structure 
of such states, we need to use formulations beyond the 
present framework, because the continuum states are not 
taken into account in the present calculations. 



IV. SUMMARY AND CONCLUDING 
REMARKS 



We have investigated the KKN system with 7=1/2 
and J p = 1/2+ as an example of multi-hadron systems 
with anti-kaons. We have performed a non-relativistic 
three-body calculation by using the effective KN inter- 
actions proposed by Hyodo-Weise and Akaishi-Yamazaki. 
With the KN potentials and no KK interaction, the 
present calculation suggests that a weakly bound KKN 
state can be formed below the KN+K threshold energy. 
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TABLE IV: Theoretical results of the KKN state calculated with the KK interactions. The range parameters of the KK 
interactions are taken to be b — 0.47 fm and b = 0.66 fm for the HW-HNJH and the AY potentials, respectively. Energies and 
root-mean-square (r.m.s.) radii and distances are listed. The expectation values for the various isospin configurations, and those 
for the unsymmetrized wave function $ are also shown. The detailed descriptions are given in captions of Tables ITT1 and HITl In 
the calculations for HW-HNJH((5tj = MeV) with U' k = k = 70 MeV, we take the extended basis space as (n max , A max ) = (15, 28) 
and (r m in, r max , -Rmin, -Rmax) = (0.2, 20, 0.2, 400) (fm) for the configuration c = 1 to get the convergent solution. 







AY+ 


Vrk 


HW-HNJH + Vkk 










5u> 


= 


5oj = 


= 11 


O-KK ( fm ) 




1.0 


1.4 


1.0 


1.4 


1.0 


1.4 


(MeV) 




70 


104 


205 


313 


205 


313 




KKN (I 


= 1/2) state 










-B(KKN) (MeV) 




-33.2 


-32.3 - 


-10.5 


-10.4 


-11.5 


-11.4 


—B(KKN) + B(KN) (MeV) 


-2.6 


-1.7 


-0.1 


-0.005 


-0.1 


-0.01 


E Im (MeV) 




-25.1 


-23.8 - 


-21.7 


-21.0 


-22.8 


-22.0 


<ImV /=0 ) (MeV) 




-22.9 


-22.1 - 


-21.5 


-21.0 


-22.5 


-22.0 


(ImV /=1 ) (MeV) 




-2.2 


-1.7 


-0.2 


0.0 


-0.3 


-0.1 


Tr (fm) 




2.7 


3.2 


12 


49 


10 


32 


expectation values for ^ 


li{[KK] ) 




0.13 


0.14 


0.20 


0.24 


0.20 


0.23 


a([KK]i) 




0.87 


0.86 


0.80 


0.76 


0.80 


0.77 


rit\iKK=o (fm) 




3.4 


3.9 


13 


50 


11 


33 


r n\iKK=i (M 




2.6 


3.0 


12 


48 


9.8 


31 


dkR (fm) 




4.1 


4.9 


18 


76 


16 


49 


expectation values for $ 


n([K 2 N] )m 




0.85 


0.87 


0.93 


0.99 


0.92 


0.98 






0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


da 2 N (fm) 




1.2 


1.3 


1.9 


2.0 


1.8 


1.9 


d (K 2 N)-R 1 (fm) 




4.2 


5.0 


19 


76 


16 


49 



The AY potential for the KN interactions provides a 
deeper bound state of the KKN system because of the 
stronger attraction in the KN channel with 1 = than 
the case of the HW potential. 

Investigating the wave function obtained by the three- 
body calculation, we have found that the KKN bound 
state can be interpreted as a hybrid of two configura- 
tions; In the inner region, two kaons are spatially sym- 
metric and couples to Ikk — 1, while in the asymptotic 
region where one kaon far from the nucleon, the system 
is regarded as the A(1405)+-A like cluster. The root- 
mean-square radius of K distribution is found to be a 
large value due to the A(1405)+AT like component with a 
loosely bound kaon around the A(1405). This kaon-halo 
like structure is reminiscent of the neutron-halo observed 
in unstable nuclei. 

We have also evaluated the decay width of the KKN 
system to irYK, obtaining T = 40 ~ 60 MeV, which 
is comparable to the width of A(1405). The dominant 
mode is found to be KKN — > K(wYl)j = o, reflecting the 
A'+A(1405) cluster structure. In this estimation, we have 
not considered three-body forces nor transitions to two- 
hadron decays. Nevertheless, such effects are expected to 
be suppressed because the overlap of the wave functions 
of three particles in a compact region is very small in 



the present case that one of the Ks is loosely bound and 
distributes very widely around the KN subsystem. 

To estimate the effect of unknown KK interactions, we 
have introduced a repulsive interaction with 1=1 which 
reproduces the K + K + scattering length obtained by lat- 
tice QCD calculation. It is interesting that a bound state 
of the KKN system is possible even with the repulsive 
KK interactions. The repulsive nature of the KK in- 
teraction suggests that the A(1405)-|-A > cluster develops 
more and the anti-kaon is further loosely bound. 

The peculiar structure that the wave function of one of 
the anti-kaon spreads for long distance is a consequence 
that the KN interaction is strong enough to form a qua- 
sibound KN and KKN states, but is not so strong for 
the deeply bound KKN state. This is very significant 
for the multi-kaon system. In such systems, some of the 
anti-kaon could be bound very loosely, so that many- 
body absorptions of anti-kaons are suppressed and the 
anti-kaons keep their identity. 

So far the experimental data for the S — — 2 channel 
have been very poor. In near future, detailed experimen- 
tal investigations could clarify the structure of excited 
5 baryons, for instance, in double strangeness reactions 
at J-PARC. The KKN molecular state suggested in the 
present calculation is one of the excited 5 baryon with 



10 




( ;- 2fm) 6.2 fin J 
(3.2 fm) ; 

FIG. 3: Schematic figure for the KKN state. K\(r) and 
K2{r) indicate the anti-kaon wave functions. One of the anti- 
kaons distributes in a wide region and the other anti-kaon 
distributes near the nucleon. The values shown in the figure 
without the parenthesis are an d d(K 2 N)-K 1 calculated 

with the HW-HNJH potential (5u> = 11 MeV), while the val- 
ues in the parentheses are those for the AY potential. In the 
inner region (I), two anti-kaons are coupled to isospin sym- 
metric with the isospin configuration [[KK]i KK = i N\ I= i/2 be- 
cause they occupy the same orbit and are spatial symmetric. 
In the outer region (II), the nucleon and one of the anti-kaons 
form a A(1405) state and the other anti-kaon is moving around 
the A(1405). 

J p = 1/2+ sitting around 1.9 GeV and having charac- 
teristic properties. The main decay mode is the three- 
body decay of the KitY, with I = for the final 7rE, 
since the A(1405) component in the three-body system 
is a doorway for the decay. In addition, productions of 
such a multi-hadron system may be strongly dependent 
on transfered momentum, since the loosely bound state 
has a large spatial distribution, which leads a softer form 
factor. These would be good indications for identifying 
the molecular state in experiments. 



In the present work, we have discussed the bound 
KKN state obtained by the the single-channel three- 
body calculation, where effects of coupled-channel 
meson-baryon interactions are taken into account as ef- 
fective single- channel KN interactions. Such coupled- 
channel effects could debase clear resonance shape for the 
three-body quasi-bound state in the spectrum, or could 
push the quasi-bound state up to the two-body contin- 
uum and the quasi-bound state could be a virtual state. 
Nevertheless, in principle, the spectra which will be ob- 
served in experiments have information of resonances and 
virtual states. For complete understanding of the KKN 
molecule state, further detailed studies involving dynam- 
ics of three-body resonances are necessary. For instance, 
coupled-channel calculations of the three-body system in- 
cluding KN «-> 7rS [42], such as Faddeev type calcula- 
tions done for the ppK~~ system in Refs. [ll|, Qjl] and 
for two-meson and one baryon system with S = — 1 in 
Refs. [43l ]. will be required to obtain qualitatively precise 
values of the energy and width. 

Our suggestion of the KKN quasibound state implies 
a possible existence of an excited baryon with a molecule 
structure of two mesons surrounding a baryon. The 
present investigation may be an important step to lead 
fundamental information on the physics of multi strange 
systems, such as anti-kaons in nuclear medium. 
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